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Abstract 

Some extensions of the Standard Model suggest existence of vectorlike fermions 
whose mass M can be by up to several orders of magnitude larger than the electroweak 
scale v , and with the Higgs mechanism providing only a small part m of this mass that 
is a fraction of v. It is shown that in a certain range of the parameters m and M 
there exist stable globules of degenerate fermionic matter made from the lighter of 
the neutral vectorlike leptons (superheavy neutrinos) with a number density much 
larger than rrihV 2 , with m/j being the Higgs boson mass. The stability is due to the 
deformation of the Higgs held inside the globule to values greatly exceeding v. The 
globules are absolutely stable if the lighter of the vectorlike neutrinos is stable and can 
be (a part of) the dark matter. 


An addition of new heavy vectorlikc families of quarks and/or leptons to the Standard 
Model generally does not lead to contradictions with the existing constraints from exper¬ 
iment. The vectorlike fermions can have a mass term that is a singlet under the elec- 
troweak (EW) gauge group 577(2). Thus such fermions are allowed to acquire ‘a mass 
of their own’, not constrained by the condition of absence of a strong Yukawa coupling 
to the Higgs boson, so that their mass scale M can be much larger than the EW scale 
v — (Gf y/2)^ 1 ^ 2 ^ 250 GeV. A possible weak Yukawa interaction with the Higgs held pro¬ 
vides a relatively small splitting m between the masses of these fermions. (In particular, 
increasing the Yukawa coupling above values corresponding to m ~ 100 GeV can create 
problems with the Higgs vacuum stability [I].) Models containing vectorlike fermions were 
considered in the literature on numerous occasions either in the context of purely theoretical 
constructions such as a unification of gauge couplings 121 Eli, or in connection with their 
specific phenomenological signature: the muon g — 2 anomaly ME], the precision elec- 
troweak measurements |8j and in the Higgs boson decays mmmmmmmmm- 
A comprehensive survey of models with vectorlike fermions can be found in Ref. [18] and a 
recent analysis of limits on such models from LHC data is presented in Ref. [19]. It can be 
also noted that if vectorlike fermions carry a conserved global quantum number the lightest 
of them, most likely electrically neutral superheavy neutrinos, constitute a part or all of the 
dark matter. 

If vectorlike superheavy neutrinos do indeed exist and were at some point present in a 
significant abundance in the Universe, it makes sense to consider their collective behavior 
at a finite density. In what follows the energy of the ground state with a uniform number 
density n of the lightest seperheavy neutrinos is considered in a model with vectorlikc leptons. 
As a result it is found that in a significant part of the parameter space of such models the 
neutrinos may form globules with the density that at large M scales as n oc / M 3 m 2 /v 2 with 
a small factor / being a function of the ratio m/rrih , where nih fs 125 GeV is the mass of 
the Higgs boson. The stability of such object with respect to dissociation into individual 
fermions arises from a deformation of the Higgs held h by the Yukawa interaction with the 
superheavy neutrinos to a value greatly exceeding its vacuum average: in the limit of large 
M the Higgs held inside the globule scales as h oc g Mm 1 / 3 /v 1 ^ 3 with a numerically small 
factor g that depends on m/m^. The enhancement of the Higgs held inside the condensate 
lowers the effective mass of the neutrino and provides a reduction of the overall mass of 
the globule below the sum of the masses of the constituent fermions in vacuum. A detailed 
consideration of possible dynamics of the process of formation of the globules is beyond the 
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scope of the present study. It will just be argued that even if a sufficient density of the 
neutrinos is created, there is still an energy barrier in the deformation of the Higgs held from 
its vacuum value v. Thus the formation of the discussed condensed objects in the presently 
existing vacuum is very improbable. However, if, in the course of cosmological evolution, 
the Higgs held in the vacuum had evolved from larger values, the globules could well be left 
over until the present epoch. 

For concreteness the treatment here is limited to a minimal model with one generation 
of vectorlike leptons with a conserved fermion number. The model thus contains two EW 
SU( 2) doublets and at least two singlets which can be labeled as being of opposite chirality: 
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with N\ and N 2 being electrically neutral and E - charged. Naturally, charged singlets can 
be added to the model, but the minimal scheme with only neutral singlets is sufficient for 
the purpose of the present paper. Moreover, a further restriction of exact L -H- R symmetry 
(the parity conservation in the vector like sector) is assumed in the present treatment, so that 
the mass term and the Yukawa interaction with the Higgs doublet (j) is written as 


A £ = —Mi (LR) - M 2 (N 2L N 2R ) — y (f> (LN 2R ) + (RN 2L ) + h.c. 


( 2 ) 


where Mi and M 2 are the EW singlet mass parameters, and y is the Yukawa coupling. In 
this scheme the charged fermion acquires mass Mi, while the mass term for the two Dirac 
neutral fermions N\ and N 2 has the form 


C N = -Mi {NiNi) 


M 2 (N 2 N 2 ) -m( 1 + x) (WW) + (N 2 N! 


( 3 ) 


where m = y/2 yv, and the notation x = h/v is introduced in a constant Higgs held, whose 
neutral component is parametrized in terms of the vacuum mean value v and the deviation 
h as 0° = (v + h)/y/2. Clearly, the mass eigenvalues for the two neutrinos are M± = 
M ± ijA 2 + m 2 ( 1 + x) 2 where M = (Mi + M 2 )/2 and A = (Mi — M 2 )/ 2, and the lighter 
of the two neutrinos is stable with respect to weak decay within this model. It is for this 
lightest of the two neutrinos that a state with a finite number density n is considered here. 

Let us start with considering the limit A = 0, so that the mass of the lighter neutrino 
in the Higgs vacuum is Mn = M — m, and later discuss the case of nonzero A. In a 
degenerate fermion state with number density n the neutrinos fill up the momentum space 
up to the Fermi momentum p, such that p 3 /(3n 2 ) = n. Thus the particles are nonrelativistic 
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inasmuch as n C M 3 , which condition, as will be seen further, is satisfied in all the range 
of M and m, where the discussed effect takes place. Allowing for the deviation h of the 
Higgs held inside the medium to be nonzero, one can write the expression for the difference 
8e(n, h) = e(n, h ) — [.M — m(h = 0 )]n between the energy density of the medium and the 
sum of the masses of the neutrinos at the density n: 


Se(n, h ) 


m 2 h v 2 2 / y \ 2 (3tt 2 n ) 5 / 3 n 2 

2 X V 2 J 10?r 2 M 4u 2 (l + y ) 2 
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(4) 


The discussed effect of the neutrino condensation is driven by the last term in this formula 
describing the shift in the neutrino mass under the shift of the Higgs held h = yu which is 
the only one term that can be negative. The hrst term on the right hand side describes the 
energy density of the shifted Higgs held, and the second term is the density of the kinetic 
energy of the degenerate nonrelativistic neutrinos, where it is assumed that m( 1 + y) M. 
Finally, the quadratic in the density third term on the right hand side in Eq.(j4]) arises from 
the self-interaction due to the Z boson force. Indeed, a matter with a finite density of 
neutrinos has a density of ‘charge’, proportional to n, with respect to the Z boson, and 
the third term describes the positive energy due to the effective Fermi self interaction of a 
neutral current, where the shift of the Higgs field mean value is also taken into account. Also 
the coefficient in this term refiects the fact that at A = 0 the lighter neutrino is an equal 
mixture of an SU( 2 ) singlet and a doublet, so that the density of the Z charge is generated 
by only one half of the number density of the particles. 

The condensation of neutrinos into globules is possible if 5e is negative at some values 
of n and y, in other terms if the negative last term in Eq.((4j) overcomes the sum of the 
three positive terms. One of the latter terms, the one with the kinetic energy, can be made 
arbitrarily small by (temporarily) considering the limit M —> oo. In this limit the condition 
for 8e to go negative does not depend on M, and it is a straightforward elementary algebra 
exercise to establish that the minimum of 5e is indeed negative if the mass parameter m is 
sufficiently large: 

m > pa 44 GeV . (5) 

2y2 

Provided that this condition is satisfied, the minimum of 5e stays negative at a sufficiently 
large, but finite, mass M. In order to estimate the asymptotic behavior of the density and 
the Higgs field at large M one can notice that the kinetic energy term does not depend on y, 
so that a minimization of the energy density with respect to y at a fixed n does not depend 
on M. At n m^v 2 one readily finds that at the minimum the value of y is also much 
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larger than one and is given by 
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with the definition r = \Jl + m 2 /m\ + m/rrih- Upon substitution in Eq. (141) this gives the 
dependence of the minimal energy density on n and M: 
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One can readily notice that the condition (J5J) is exactly the one for the coefficient of the 
n 4 / 3 term to be negative. Once this condition is satisfied the expression in Eq.([7|) develops 
a negative minimum at the density n such that 


n (8rm/rrih—2 — r 2 ) 3 rrihM 3 
rrihV 2 243 7r 4 r 2 v A 


( 8 ) 


Clearly, the proportionality of n to M 3 justifies at large M the assumption n rrihV 2 , used 
in the derivation. On the other hand, Eq. (jHJ) gives n M 3 , which justifies the nonrelativistic 
limit used in the evaluation of the kinetic energy of the neutrinos. Also, according to Eq. (j6jl . 
the deformation of the Higgs mean value inside the condensate is strong, y 1, but still 
X -C M/m, so that the shift of the neutrino mass by the Yukawa interaction is small as 
compared to M. 

Having established the effect of the condensation of the (superheavy) neutrinos at a 
large M, it is natural to address the question of the minimal value of M at a given m , 
satisfying the condition (]5]h at which the condensation is possible. Near the critical value of 
M the value of y at the minimum of the energy density is generally not very large, the finite 
terms in the y dependence in Eq. (HI) can not be neglected, and the minimization is better 
done numerically. In performing such minimization it is helpful to first minimize the energy 
density with respect to n. Indeed, the expression in Eq.(j4l) is polynomial in n 1//3 and the 
minimization in n is equivalent to solving a cubic equation, which can be done analytically. 
The remaining after that minimization in y can be readily done numerically. The result of 
such numerical evaluation is illustrated in Fig. 1, where the minimal value of M, above which 
the condensation takes place, is shown as a function of m. As it follows from the previous 
discussion the minimal value of M approaches infinity when m goes down to the value in 
the r.h.s. of Eq.(l5il. 

It should be also noted that at any given n the energy density 5e is always positive at 
y -C 1 which can be readily verified by considering the expansion of the expression in Eq. (J4j) 
to the second order in y. Thus a formation of the discussed globules always encounters a 
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Figure 1: The dependence on the parameter m of the minimal value of the mass M, above 
which the condensation of heavy vectorlike neutrinos into globules is possible. 


barrier before the Higgs held mean value is deformed by at least h ~ O(v). The scale for 
the height of barrier is set by Higgs energy density ~ m|n 2 . Therefore in any reasonable 
situation at present the formation of the globules is essentially impossible. 

The expression in Eq.(j4j) for 5e assumes, as previously mentioned, that the mass term 
(l3j) contains the degeneracy Mi = Mo, so that A = 0. For a nonzero A one can consider two 
cases: when A is comparable to M Q], and the case A M. The first case in fact splits into 
two sub-cases corresponding to different sign of A. If A < 0, the lighter neutrino is mostly a 
component of an SU(2 ) doublet and the positive contribution to the energy density due to 
the Z boson held enters with a four times larger coefficient than in in Eq. (J4J). On the other 
hand, as long as A ^ m(l + y), the effect of the mass dependence on the Higgs mean held is 
approximately equivalent to replacing in Eq. (J4]) the parameter m by m 2 (l + y/2)/A <C m, 
so that the formation of the condensate is impossible since the effective parameter would 
violate analog of the condition ((5]). If however m(l + y) ~ A ~ M, the neutrinos inside the 


lr The case A = M is the well known l see-saw’ mechanism for the masses of the ordinary (chiral) neutrinos. 
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globule have to be relativistic, in which situation the effect of the (weak) dependence of mass 
on the Higgs field mean value cannot compensate the positive kinetic energy of the particles. 
The sub-case where A > 0 is potentially more interesting, since the lighter neutrino is then 
mostly a singlet, and the self-repulsion through the Z boson is greatly reduced. In the limit 
m(l + x) A the difference Se of the energy density takes the form 
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It can be readily verified that the necessary condition for this expression to develop a negative 
minimum is m > y/2mh ~ 177 GeV, which would be problematic from the point of vacuum 
stability [lj. 

Thus only the case A <C M appears to be interesting with regards to the existence of 
the condensate. In this case for a sufficiently small A one can find m(l + x) > A so that 
the effect of the mass splitting A is small. Using the formulas in Eq.flB]) and Eq.flH]) it can 
be estimated that this regime is achieved if A/M is much smaller than ( m,h/v ) 4//3 , if m and 
rrih are considered as being of the same order. 

The presented estimates are based on considering the energy density of a spatially uniform 
condensate of the heavy neutrinos. Such treatment is justified only inasmuch as the surface 
energy can be neglected in comparison with the bulk energy, associated with the transition 
of the Higgs field on the surface of the globule, and requires the radius of the globule R to be 
large. Considering the parameters m and rriy, as being of the same order and describing the 
smallest mass scale in the problem, one can conservatively estimate that the condition for 
the approximation used is R rrih 1. According to the estimate in Eq. (JHJ) for the equilibrium 
number density at large M, the total number AT of particles in the globule should then satisfy 
the condition 

4(8 r m/rrih — 2 — r 2 ) 3 


AT » 
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( 10 ) 


729 7r 3 r 2 nth v 2 

The estimate (TTOT) of the total number of particles naturally places a limit on how big the 
mass M can be before the globule collapses into a Black Hole. A straightforward estimate 
then gives the condition 

M < y/Mpiv ~ 10 11 GeV , (11) 

still leaving an ample range of the values of M at which the discussed globules may exist. It 
is also interesting to notice that within this range the scale for the mass MJ\T of the globule 
can reach macroscopic values: 

M Y 
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while the scale for its radius is set by 1/rrih ~ 10~ 16 cm. 

The possibility that the discussed here globules do indeed exist makes it very interesting 
to study their phenomenological implications, in particular their interaction with ordinary 
matter. Any detailed discussion of this subject is beyond the scope of the present paper. 
However, it can be mentioned, as a simple remark, that the enhanced Higgs held inside the 
globule results in a strong repulsion from the ordinary matter. Indeed, the masses of leptons 
and quarks are increased inside the condensate by the factor (1 + x) and the QCD mass scale 
Aqcd is larger inside the globule due to the larger masses of quarks and the resulting cutoff 
at a higher scale of the screening contribution to the running of the QCD coupling. For this 
reason it is impossible to exclude existence of the globules from the apparent absence of such 
objects bound in ordinary matter. 

In summary. In an extension of the Standard Model with vectorlike leptons it is found 
that in a certain part of the parameter space the lighter (but still superheavy) neutral 
fermions can assemble into globules with the number density n exceeding the EW scale 
mh,v 2 . The energy per particle inside the globule is lower than the fermion mass in the 
vacuum due to the strong enhancement of the Higgs field inside the condensate by the 
collective effect of the Yukawa interaction. The mass of a globule can reach a macroscopic 
scale (i.e. up to grams or more), and the globules are absolutely stable if the constituent 
neutral fermions are stable and can be a part or all of the dark matter. However, if at present 
the vectorlike fermions exist in a form of dispersed individual particles, their spontaneous 
condensation into globules is prevented by an energy barrier. If the globules could still be 
produced during the cosmological evolution, a study of their properties, such as interaction 
with ordinary matter, can result in new limits on models with vectorlike fermions, or in a 
realization of a possible new form of the dark matter. 
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